gamma=0.1;
r=0;
gammaBase=gamma;
rBase=r;
face=1;

Lambda=1;
mu=1;
kappa=9;
phi=9;


tolerance=0.0001;
Agency=1;
minimal=0;
delta=0;



aBar=0.125;
qBar=0.24;
cond=phi*aBar/(gamma+Lambda-qBar-aBar)/qBar;
cond0=2*mu/(gamma+Lambda-qBar-aBar)^2/(r+Lambda-qBar-aBar);
Vmax=1/(gamma+Lambda-qBar-aBar)/(r+Lambda-qBar-aBar);
cond2=(2*Vmax+Vmax^2/phi)/(gamma+Lambda-qBar-aBar);


qGrid=200;

options = odeset('AbsTol',1e-7,'RelTol',1e-7);


alpha=0;
alphaNeg1=-2;
alphaNeg2=-1;


K=0.75;
pLambda=0;
R=0;